Acoustic measurement of the low-energy excitations in Nd2-xCea;Cu04+5 
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The complex dynamic Young's modulus of ceramic Nd2-iCei,Cu04 with s = 0, 0.05 and 0.20 
has been measured from 1.5 to 100 K at frequencies of 1 — 10 kHz. In the undoped sample the 
modulus starts decreasing below ~ 20 K, instead of approaching a constant value as in a normal 
solid. The modulus minimum has been interpreted in terms of paraelastic contribution from the 
relaxation of the Nd^+ 4/ electrons between the levels of the ground state doublet, which is split 
by the interaction with the antiferromagnetically ordered Cu sublattice. The value of the splitting 
is found to be 0.34 meV, in excellent agreement with inelastic neutron scattering, infrared and 
specific heat experiments. With doping, the anomaly shifts to lower temperature and decreases in 
amplitude, consistently with a reduction of the local field from the Cu sublattice. 
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I. INTRODUCTION 

The crystal-field (CF) excitations of the Nd^+ 4/ elec- 
trons in Nd2-a:Cea;Cu04 have been extensivejv-studied, 
mainly by inelastic neutron scattering (INS)Eru and in- 
frared spectroscopies^ Although discrepancies exist be- 
tween the proposed sets of phenomenological parameters 
in the CF hamiltjonian, there is general agreement on the 
level scheme.ElD'El The ground state multiplet of Nd'^+ 
has J — 9/2 and is split into five doublets, the sepa- 
ration of the first two being 14 meV; each doublet, in 
turn, is split by the interaction of the Nd magnetic mo- 
ment with the local magnetic field, mainly created by the 
antiferromagnetic ordering of the Cu magnetic moments 
below T]^" ~ 260 K. The magnetic separation of the first 
doublet has been determined by INS, and passes from 
0.39 meV for the undoped case to 0.2 meV for x — O.IS.eI 
Such a low-energy excitation also appears as a|-Seliottky 
anomaly in the low-temperature specific heat,Enl3 with 
values of the level splitting in agreement with the INS re- 
sults; in principle, it should also give rise to an anomaly 
in the elastic moduli, if the level separation can be modu- 
lated by strain. To our knowledge, only two experiments 
have been reported on the-elastic properties of NCCO at 
liquid He temperatures Jlilll3 but no attempts have been 
made to relate those data to the elastic response of the 
ground state doublet. We made anelastic spectroscopy 
measurements in the kHz range on Nd2-2:Cea;Cu04 with 
a; = 0, 0.05 and 0.20, and found a decrease of the Young's 
modulus at liquid He temperatures, which can be iden- 
tified with the paraelastic response of the ground state 
doublet of Nd3+. 



II. EXPERIMENTAL 

The samples were prepared by solid state reaction, 
with repeated sintering and grinding of the powders up 
to obtaining X-ray diffraction spectra free of unwanted 
peaks; the methods used, together with a thorough ther- 
modynamic study of the phase diagram and stability of 
the Nd-Ce-Cu-0 system, are reported in Ref. Gen- 
erally, a reducing treatment is made on Nd2-2;Cea;Cu04, 
in order to make it superconducting; the data presented 
here refer to the as-prepared or reoxygenated status. 

The samples were cut as bars approximately 40 x 4 x 
0.6 mm^. The dynamic Young's modulus, E{lu,T) — 
E' + iE" , was measured by suspending the bars on their 
nodal lines on thin wires and electrostatically exciting 
their flexural vibrations; the 1st and 5th flexural modes 
were chosen, whose resonance frequencies iojlix are ap- 
proximately in the ratio 1 : 13.3. The vibration am- 
plitude was measured by a frequency modulation tech- 
nique, where the capacity between the sample and the 
electrode is inserted into a high frequency oscillator, and 
the resulting signal, after passing through an amplitude- 
modulation discriminator, is sent to a lock-in amplifier. 



The real part E' of the Young's modulus is given by 
E = p(^0.9752^^^ , where p, I and h are the sample 

density, length and thickness; a correction for the poros- 
ity should be necessary, but we are only interested in the 
relative variation of E with temperature. The elastic en- 
ergy loss coefficient = E" /E' was measured from 
the free decay of the vibration amplitude after switching 
off the excitation signal. 



III. RESULTS 

Figure 1 presents the complex Young's modulus of 
Nd2-a;Ce£i;Cu04 at three doping levels, between 1 and 
100 K; only the data of the first mode are presented, since 
those at higher frequency are identical. The elastic en- 
ergy loss coefficient decreases monotonously with temper- 
ature and with increasing doping, without any particular 
feature except for a small peak at 15 K for the undoped 
sample. Such a peak disappears after annealing in vac- 
uum at high temperature, and, although several reorien- 
tational transitions of the Nd and Cu spins have beejMb-. 
served at low temperature with various techniques 
it does not seem to be associated to any of these mag- 
netic transitions. We will not discuss this peak further. 
The elastic modulus presents a nearly linear stiffening on 
cooling below room temperature, which levels off below 
30 K; contrary to a normal solid, however, the modu- 
lus does not stabilize on further cooling, but drops of 
several parts in 10'^. The drop actually saturates at the 
lowest temperatures reached in the experiments (see also 
Fig. 2), and in some cases we measured the beginning 
of a rise below 1.5 K. The temperature and amplitude of 
the drop decrease with increasing doping similarly to the 
Schottky anomaly in the specific heat, and therefore this 
feature might well be connected with the ground state 
doublet of the Nd'^^ ions. 



IV. DISCUSSION 

Most of the specific heat data below 10 K have been 
interpreted in terms of a Scbottky anomaly arising from 
the ground state doublet J£ril3 with a splitting A (T) = 
2//o^Cu (T) determined through Nd-Cu exchange by the 
staggered magnetization Bqu {T) of the AF Cu sublattice 
below Tj^^"; since Tj^" ~ 260 K, Bcu and the splitting A 
are practically constant below 50 K. The existence of two 
levels Ei separated by A, whose internal energy is U = 
(A/2) tanh (A/2A:Br), contributes to the molar specific 
heat with the Schottky term 



SCj, ~ — — — knc ( — — — I ^ 

" dT ^ \2kBT) cosh^ (A/2fcBT) 



(1) 



where c is the molar concentration of the relaxing entities, 
in the present case the concentration 1 — a;/2 of Nd ions. 
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On the other hand, if the sphtting cam jje modulated by 
a strain e, paraelastic and diaelastictSEZl contributions to 
the elastic modulus M = d'^U/de^ are expected. Here 
we omit the tensoriai character of the elastic proper- 
ties, since we are dealing with the Young's modulus of 
polycrystalline samples, which contains a combination of 
several elastic constants and we are not entering in the 
details of the effect of the strain symmetry on the Nd'^^ 
multiplets; then M is identified with the Young's modu- 
lus E and e is the corresponding effective strain. i— ii— i 
The paraelastic contribution can be written aalj'tZl 
SMpara = z2ij If" ^ "aT *° relaxation 

of the populations n;, perturbed by the vibration strain 
e through the modulation of the eujecgies Ei. For only 
two levels E2 — Ei — A, one obtainsO 
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4fcBrcosh2 (A/2fcBr) ' 



(2) 



where vq is the molecular volume and 7 = dA/de is 
the deformation potential. Such a relaxation occurs 
with a cba^acteristic rate t^^, and therefore the dynamic 
modulusta has to be multiplied by a factor (1 -I- zwr)~^. 
The resulting real part of the dynamic modulus M' (lo, T) 
acquires a frequency-dependent dispersion, while the 
imaginary part produces absorption in correspondence 
with the modulus dispersion. The present data, however, 
do not exhibit any frequency dependence of M' (w, T) or 
absorption peak, and this implies ojt ~ 0, namely the 
rate is much faster than the measuring frequencies 
uj ~ 10^ — 10^ s""^ . Indeed, the Cu and Nd spiiwates de- 
duced from muon spin relaxation experiments^ in sam- 
ples with X = 0.20 are faster than 10^" s^^. 

If the first derivative 7 of the energy split with respect 
to strain is null for symmetry reasons, than one remains 
with the diaelastic term, cJM'^''' = — y,,n^^^^, which 
contains the second derivatives. The physical meaning of 
such a term is simply that the curvature of each energy 
level versus strain contributes to the corresponding elas- 
tic constant, and its contribution is weighted with the 
level population. For only two levels one obtains:t3 



cld^A , 
<J Afdia = — - „ o tanh 



vq 2 de 



2kBT 



(3) 



The continuous curves in Fig. 2 are fits with (SAfpara 
given by Eq. |^ and assuming that the background mod- 
ulus depends on temperature like an insulating Debye 
solid,E3 M{T) = Mo(l-aT4). An additional term 
quadratic in T is expected when the system becomes 
metallic, due to the energy of the free electrons, which 
increases quadratically with T; such a term, however, 
could be of importance only for the highest doping 
X = 0.20, and still at x = 0.15 some elastic constants 
have been found to ba-well described by the quartic 
term alone below 60 K.t3 From the data at cc = and 
0.05 we find that M (T) is described by the coefficient 
a = -2.1 X lO"^'^ K"''. The values of the splitting de- 
duced from the fits of Fig. 2 are A/k-Q ~ 4.0, 2.25 and 



1.3 K at x = 0, 0.05 and 0.20 respectively, while the 
amplitudes of the paraelastic contributions are in the ra- 
tios 1 : 0.6 : 0.05. The sphtting found for the undoped 
case is inj-jexcellent agreement with the values deduced 
from INSB (A/fce = 41)6 K ), infrared spectroscopy^ 
(3.8 K) and specific heatQ (A/Zcb = 4.5 K). Doping by Ce 
substitution introduces electrons mainly of Cu 3d char- 
acter, which cancel the magnetic moments of those Cu 
ions and therefore reduce the CitNd exchange responsi- 
ble for the ground state splittingjj As a result, the Schot- 
tky anomaly in the specific heat has been found to shift 
to lower temperatua;_|tp broaden and to decrease in in- 
tensity with doping.Q'Ba The present data confirm such a 
trend, but the reduction in intensity is more marked than 
in the specific heat measurements possibly because of a 
reduction of the strain dependence of the splitting with 
doping. 

At x > 0.15 there is no long range AF order of the 
Cu spins, and the magnetic ground state splitting should 
disappear; nonethjcless, the specific heat presents a peak 
up to X = 0.20Jij Our data &t x — 0.20 also present a 
residue due to the doublet with reduced splitting, but 
the value deduced from the fit, A/k-Q — 1.2 K, should 
be considered with caution in view of the possible influ- 
ence of an electronic quadratic term in the temperature 
dependence of the modulus. 

In some cases the fits to the specific heat data have 
been improved by assuming an additional temperature 
dependence of the magnetic, splitting attributed to the 
ordering of the Nd momentaJj'Ll which however affects the 
splitting by less than 10%, or by assuming a three-level 
structure instead of a doublet. El We chose to not over- 
interpret our data with additional parameters besides a 
simple Schottky anomaly with constant splitting below 
50 K and the dependence of the background elastic 
modulus. In fact, the contribution of the doublet to the 
elastic response, Eq. ^ has a slower temperature depen- 
dence than the contribution to the specific heat, Eq. (j^), 
in the high temperature limit, being proportional to 
instead of T~^; therefore, the elastic anomaly is spread 
to higher temperature where additional effects become 
important. In addition, the lowest temperature attained 
in our experiments was not sufficient to complete the low- 
temperature side of the curves, so that some uncertainty 
in the determination of the minimum of <5Afpara (T) ex- 
ists. This means that, while it is easy to improve the fit 
by adding a temperature dependence of the splitting, by 
taking into account the contribution of the excited dou- 
blets, and adopting more sophisticated expressions for 
M (T) , it is difficult to distinguish which of the improve- 
ments is the important one. Still, the data presented 
here, together with other results on reduced samples, are 
sufficient to indicate the presence of a minimum in the 
modulus due to a paraelastic response, rather than a 
step-like anomaly due to a diaelastic response. Fits to 
the dialeastic expression, Eq. |^, are less satisfactory and 
yield values of A about twice larger than those of Fig. 
2, in contrast with the literature values. If both the di- 
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aelastic and paraelastic terms are considered, with both 
7 and 9^A/i9e^ as free parameters, the splitting A re- 
mains almost the same as reported here and up to half of 
the anomaly can be accounted for by the diaelastic term, 
without a significant improvement of the fit. We conclude 
that the anelastic response in the modulus is mainly due 
to the linear part of the response of the ground-state 
splitting to strain. The magnitude of the deformation 
potential 7 can be extracted inserting in Eq. ^ the value 
obtained from the fit, {SM^^^^/M)^^^ = 0.0044 at 2 K for 
X — 0, and using A/ ~ 100 GPa, as determined from the 
resonance frequency (without corrections for the poros- 
ity); it turns out 7 = 0.013 eV for x — Q and 0.01 eV for 
X = 0.05. 

As noted above, the absence of a rise of the absorp- 
tion in correspondence to the decrease of the modulus 
implies that the relaxation rate for the Nd'^+ spins to 
redistribute themselves within the ground states dou- 
blet is much faster than our measuring frequency, lo ^ 
10-* - 10^ s-^. Also, there are no signs of freezing of the 
Cu spins into a spin-glass or cluster glass states, which 
instead is observed in the hole-doped Cu02 planes of 
La2-a;Sr2;Cu04 and YBa2Cu306+a;;" in La2-a:Sra:Cu04 
the frozen spin donpains produce a distinct rise of the 
acoustic absorption,c3 which is totally absent here. Sim- 
ilar conclusioHS have been drawn from /iSR relaxation 
experiments.llj 1— in 

To our knowledge, only two other experimental3'Lll ex- 
ist on the acoustic properties of Nd2-£!;Cea;Cu04 at liq- 
uid He temperatures, which however cannot be com- 
pared with the present one in a straightforward way; in 
fact, they report some elastic constants of single crystals, 
while we measured the Young's modulus of polycryatals, 
whickxiontain a combination of all the compliances. E£l Fill 
et alxA found several anomalies in the elastic constants 
of undoped Nd2Cu04 as a function of temperature and 
magnetic field. In particular, the cge shear presented 
a minimum at 5 K with an amplitude nearly ten times 
larger than our minimum, and with a narrower shape, 
followed by a smaller decrease below 1 K tentatively at- 
tributed to a Schottky anomaly or magnetic ordering of 
Nd ions. Around 5 K also step-like changes have been 
found in other elastic constants, which have been inter- 
preted as a fepcnalastic transition driven by an ordering of 
the Nd spins O'Ej We cannot find a clear correspondence 
between those data and the ones presented here. |— ■ 

Regarding the doped material, Saint-Paul et a/.EJ 
measured the sound velocity of a single crystal of 
Ndi.85Ceo.i5Cu04; the shear mode in the ab plane did 
not present any anomaly down to 10 K, but this is not in 
contrast with our data at x = 0.20, where the decrease of 
the modulus is hardly detectable above 10 K. Instead, the 
C33 mode exhibited an upward deviation below 10 K with 
respect to the —T'^ dependence, opposite to our measure- 
ments. A possible explanation for the discrepancy is that 
the major contribution to the C33 mode is diaelastic with 
positive 9^A/9e^, while the Young's modulus of the poly- 
crystal contains a predominant paraelastic contribution 



which cancels the diaelastic one. 



V. CONCLUSION 

The Young's modulus and elastic energy loss coeffi- 
cient of Nd2-a;Cea;Cu04 with x = 0, 0.05 and 0.20 have 
been measured from 1.5 to 100 K. On approaching the 
lowest temperature, the modulus presents a drop whose 
temperature and amplitude decrease with increasing dop- 
ing. The data have been interpreted in terms of parae- 
lastic contribution from the relaxation of the Nd'^"'" 4/ 
electrons between the levels of the ground state dou- 
blet, which is split by the interaction with the antifcr- 
romagnetically ordered Cu^^ ions. Excellent agreement 
is found with the value of the splitting at zero doping de- 
duced from INS, infrared and specific heat measurements, 
and the shift of the anomaly to lower temperature with 
doping is consistent with a reduction of the local field 
from the Cu sublattice. An effective deformation poten- 
tial 7 = 0.013 eV is found for the strain derivative of 
the splitting at x = 0, which decreases with doping. The 
absorption is very low and featureless, indicating that 
the relaxation rate for reaching equilibrium within the 
ground state doublet is much faster than the measuring 
frequency, r"^ ^ 10^ s~^, and also excluding the occur- 
rence of freezing phenomena of the Cu or Nd spins down 
to 1.5 K. 
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FIG. 1. Elastic energy loss coefficient and relative change 
of the Young's modulus E of Nd2-iCe2,Cu04 samples with 
x = Q (0.9 kHz), X = 0.05 (2 kHz) and x = 0.20 (1.6 kHz). 

FIG. 2. Fit of the relative change of the Young's modulus 
of Nd2-a;Cea;Cu04 with the paraelastic contribution from a 
doublet with splitting A, Eq. ^. 
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